Systematic mapping of TF-mediated cell fate changes by a pooled induction coupled with scRNA-seq and multi-omics approaches.

Transcriptional regulation controls cellular functions through interactions between transcription factors (TFs) and their chromosomal targets. However, understanding the fate conversion potential of multiple TFs in an inducible manner remains limited. Here, we introduce iTF-seq as a method for identifying individual TFs that can alter cell fate toward specific lineages at a single-cell level. iTF-seq enables time course monitoring of transcriptome changes, and with biotinylated individual TFs, it provides a multi-omics approach to understanding the mechanisms behind TF-mediated cell fate changes. Our iTF-seq study in mouse embryonic stem cells identified multiple TFs that trigger rapid transcriptome changes indicative of differentiation within a day of induction. Moreover, cells expressing these potent TFs often show a slower cell cycle and increased cell death. Further analysis using bioChIP-seq revealed that GCM1 and OTX2 act as pioneer factors and activators by increasing gene accessibility and activating the expression of lineage specification genes during cell fate conversion. iTF-seq has utility in both mapping cell fate conversion and understanding cell fate conversion mechanisms.

DeepSLICEM: Clustering CryoEM particles using deep image and similarity graph representations.

Finding the 3D structure of proteins and their complexes has several applications, such as developing vaccines that target viral proteins effectively. Methods such as cryogenic electron microscopy (cryo-EM) have improved in their ability to capture high-resolution images, and when applied to a purified sample containing copies of a macromolecule, they can be used to produce a high-quality snapshot of different 2D orientations of the macromolecule, which can be combined to reconstruct its 3D structure. Instead of purifying a sample so that it contains only one macromolecule, a process that can be difficult, time-consuming, and expensive, a cell sample containing multiple particles can be photographed directly and separated into its constituent particles using computational methods. Previous work, SLICEM, has separated 2D projection images of different particles into their respective groups using 2 methods, clustering a graph with edges weighted by pairwise similarities of common lines of the 2D projections. In this work, we develop DeepSLICEM, a pipeline that clusters rich representations of 2D projections, obtained by combining graphical features from a similarity graph based on common lines, with additional image features extracted from a convolutional neural network. DeepSLICEM explores 6 pretrained convolutional neural networks and one supervised Siamese CNN for image representation, 10 pretrained deep graph neural networks for similarity graph node representations, and 4 methods for clustering, along with 8 methods for directly clustering the similarity graph. On 6 synthetic and experimental datasets, the DeepSLICEM pipeline finds 92 method combinations achieving better clustering accuracy than previous methods from SLICEM. Thus, in this paper, we demonstrate that deep neural networks have great potential for accurately separating mixtures of 2D projections of different macromolecules computationally.

Label-free proteomic comparison reveals ciliary and nonciliary phenotypes of IFT-A mutants.

DIFFRAC is a powerful method for systematically comparing proteome content and organization between samples in a high-throughput manner. By subjecting control and experimental protein extracts to native chromatography and quantifying the contents of each fraction using mass spectrometry, it enables the quantitative detection of alterations to protein complexes and abundances. Here, we applied DIFFRAC to investigate the consequences of genetic loss of Ift122, a subunit of the intraflagellar transport-A (IFT-A) protein complex that plays a vital role in the formation and function of cilia and flagella, on the proteome of Tetrahymena thermophila. A single DIFFRAC experiment was sufficient to detect changes in protein behavior that mirrored known effects of IFT-A loss and revealed new biology. We uncovered several novel IFT-A-regulated proteins, which we validated through live imaging in Xenopus multiciliated cells, shedding new light on both the ciliary and non-ciliary functions of IFT-A. Our findings underscore the robustness of DIFFRAC for revealing proteomic changes in response to genetic or biochemical perturbation.

Systematic profiling of ale yeast protein dynamics across fermentation and repitching.

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is among the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout 2 fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

SARS-COV-2 Omicron variants conformationally escape a rare quaternary antibody binding mode.

The ongoing evolution of SARS-CoV-2 into more easily transmissible and infectious variants has provided unprecedented insight into mutations enabling immune escape. Understanding how these mutations affect the dynamics of antibody-antigen interactions is crucial to the development of broadly protective antibodies and vaccines. Here we report the characterization of a potent neutralizing antibody (N3-1) identified from a COVID-19 patient during the first disease wave. Cryogenic electron microscopy revealed a quaternary binding mode that enables direct interactions with all three receptor-binding domains of the spike protein trimer, resulting in extraordinary avidity and potent neutralization of all major variants of concern until the emergence of Omicron. Structure-based rational design of N3-1 mutants improved binding to all Omicron variants but only partially restored neutralization of the conformationally distinct Omicron BA.1. This study provides new insights into immune evasion through changes in spike protein dynamics and highlights considerations for future conformationally biased multivalent vaccine designs.

Systematic Profiling of Ale Yeast Protein Dynamics across Fermentation and Repitching.

Studying the genetic and molecular characteristics of brewing yeast strains is crucial for understanding their domestication history and adaptations accumulated over time in fermentation environments, and for guiding optimizations to the brewing process itself. Saccharomyces cerevisiae (brewing yeast) is amongst the most profiled organisms on the planet, yet the temporal molecular changes that underlie industrial fermentation and beer brewing remain understudied. Here, we characterized the genomic makeup of a Saccharomyces cerevisiae ale yeast widely used in the production of Hefeweizen beers, and applied shotgun mass spectrometry to systematically measure the proteomic changes throughout two fermentation cycles which were separated by 14 rounds of serial repitching. The resulting brewing yeast proteomics resource includes 64,740 protein abundance measurements. We found that this strain possesses typical genetic characteristics of Saccharomyces cerevisiae ale strains and displayed progressive shifts in molecular processes during fermentation based on protein abundance changes. We observed protein abundance differences between early fermentation batches compared to those separated by 14 rounds of serial repitching. The observed abundance differences occurred mainly in proteins involved in the metabolism of ergosterol and isobutyraldehyde. Our systematic profiling serves as a starting point for deeper characterization of how the yeast proteome changes during commercial fermentations and additionally serves as a resource to guide fermentation protocols, strain handling, and engineering practices in commercial brewing and fermentation environments. Finally, we created a web interface (https://brewing-yeast-proteomics.ccbb.utexas.edu/) to serve as a valuable resource for yeast geneticists, brewers, and biochemists to provide insights into the global trends underlying commercial beer production.

Protein nonadditive expression and solubility contribute to heterosis in Arabidopsis hybrids and allotetraploids.

Hybrid vigor or heterosis has been widely applied in agriculture and extensively studied using genetic and gene expression approaches. However, the biochemical mechanism underlying heterosis remains elusive. One theory suggests that a decrease in protein aggregation may occur in hybrids due to the presence of protein variants between parental alleles, but it has not been experimentally tested. Here, we report comparative analysis of soluble and insoluble proteomes in Arabidopsis intraspecific and interspecific hybrids or allotetraploids formed between A. thaliana and A. arenosa. Both allotetraploids and intraspecific hybrids displayed nonadditive expression (unequal to the sum of the two parents) of the proteins, most of which were involved in biotic and abiotic stress responses. In the allotetraploids, homoeolog-expression bias was not observed among all proteins examined but accounted for 17-20% of the nonadditively expressed proteins, consistent with the transcriptome results. Among expression-biased homoeologs, there were more A. thaliana-biased than A. arenosa-biased homoeologs. Analysis of the insoluble and soluble proteomes revealed more soluble proteins in the hybrids than their parents but not in the allotetraploids. Most proteins in ribosomal biosynthesis and in the thylakoid lumen, membrane, and stroma were in the soluble fractions, indicating a role of protein stability in photosynthetic activities for promoting growth. Thus, nonadditive expression of stress-responsive proteins and increased solubility of photosynthetic proteins may contribute to heterosis in Arabidopsis hybrids and allotetraploids and possibly hybrid crops.

An amino acid-resolution interactome for motile cilia illuminates the structure and function of ciliopathy protein complexes.

Motile cilia are ancient, evolutionarily conserved organelles whose dysfunction underlies motile ciliopathies, a broad class of human diseases. Motile cilia contain myriad different proteins that assemble into an array of distinct machines, so understanding the interactions and functional hierarchies among them presents an important challenge. Here, we defined the protein interactome of motile axonemes using cross-linking mass spectrometry (XL/MS) in Tetrahymena thermophila. From over 19,000 XLs, we identified 4,757 unique amino acid interactions among 1,143 distinct proteins, providing both macromolecular and atomic-scale insights into diverse ciliary machines, including the Intraflagellar Transport system, axonemal dynein arms, radial spokes, the 96 nm ruler, and microtubule inner proteins, among others. Guided by this dataset, we used vertebrate multiciliated cells to reveal novel functional interactions among several poorly-defined human ciliopathy proteins. The dataset therefore provides a powerful resource for studying the basic biology of an ancient organelle and the molecular etiology of human genetic disease.

Distinctive interactomes of RNA polymerase II phosphorylation during different stages of transcription.

During eukaryotic transcription, RNA polymerase II undergoes dynamic post-translational modifications on the C-terminal domain (CTD) of the largest subunit, generating an information-rich PTM landscape that transcriptional regulators bind. The phosphorylation of Ser5 and Ser2 of CTD heptad occurs spatiotemporally with the transcriptional stages, recruiting different transcriptional regulators to Pol II. To delineate the protein interactomes at different transcriptional stages, we reconstructed phosphorylation patterns of the CTD at Ser5 and Ser2 in vitro. Our results showed that distinct protein interactomes are recruited to RNA polymerase II at different stages of transcription by the phosphorylation of Ser2 and Ser5 of the CTD heptads. In particular, we characterized calcium homeostasis endoplasmic reticulum protein (CHERP) as a regulator bound by phospho-Ser2 heptad. Pol II association with CHERP recruits an accessory splicing complex whose loss results in broad changes in alternative splicing events. Our results shed light on the PTM-coded recruitment process that coordinates transcription.

Robust and scalable single-molecule protein sequencing with fluorosequencing.

The need to accurately survey proteins and their modifications with ever higher sensitivities, particularly in clinical settings with limited samples, is spurring development of new single molecule proteomics technologies. Fluorosequencing is one such highly parallelized single molecule peptide sequencing platform, based on determining the sequence positions of select amino acid types within peptides to enable their identification and quantification from a reference database. Here, we describe substantial improvements to fluorosequencing, including identifying fluorophores compatible with the sequencing chemistry, mitigating dye-dye interactions through the use of extended polyproline linkers, and developing an end-to-end workflow for sample preparation and sequencing. We demonstrate by fluorosequencing peptides in mixtures and identifying a target neoantigen from a database of decoy MHC peptides, highlighting the potential of the technology for high sensitivity clinical applications.

Integrated modeling of the Nexin-dynein regulatory complex reveals its regulatory mechanism.

Cilia are hairlike protrusions that project from the surface of eukaryotic cells and play key roles in cell signaling and motility. Ciliary motility is regulated by the conserved nexin-dynein regulatory complex (N-DRC), which links adjacent doublet microtubules and regulates and coordinates the activity of outer doublet complexes. Despite its critical role in cilia motility, the assembly and molecular basis of the regulatory mechanism are poorly understood. Here, using cryo-electron microscopy in conjunction with biochemical cross-linking and integrative modeling, we localize 12 DRC subunits in the N-DRC structure of Tetrahymena thermophila. We also find that the CCDC96/113 complex is in close contact with the DRC9/10 in the linker region. In addition, we reveal that the N-DRC is associated with a network of coiled-coil proteins that most likely mediates N-DRC regulatory activity.

Molecular complex detection in protein interaction networks through reinforcement learning.

BACKGROUND: Proteins often assemble into higher-order complexes to perform their biological functions. Such protein-protein interactions (PPI) are often experimentally measured for pairs of proteins and summarized in a weighted PPI network, to which community detection algorithms can be applied to define the various higher-order protein complexes. Current methods include unsupervised and supervised approaches, often assuming that protein complexes manifest only as dense subgraphs. Utilizing supervised approaches, the focus is not on how to find them in a network, but only on learning which subgraphs correspond to complexes, currently solved using heuristics. However, learning to walk trajectories on a network to identify protein complexes leads naturally to a reinforcement learning (RL) approach, a strategy not extensively explored for community detection. Here, we develop and evaluate a reinforcement learning pipeline for community detection on weighted protein-protein interaction networks to detect new protein complexes. The algorithm is trained to calculate the value of different subgraphs encountered while walking on the network to reconstruct known complexes. A distributed prediction algorithm then scales the RL pipeline to search for novel protein complexes on large PPI networks. RESULTS: The reinforcement learning pipeline is applied to a human PPI network consisting of 8k proteins and 60k PPI, which results in 1,157 protein complexes. The method demonstrated competitive accuracy with improved speed compared to previous algorithms. We highlight protein complexes such as C4orf19, C18orf21, and KIAA1522 which are currently minimally characterized. Additionally, the results suggest TMC04 be a putative additional subunit of the KICSTOR complex and confirm the involvement of C15orf41 in a higher-order complex with HIRA, CDAN1, ASF1A, and by 3D structural modeling. CONCLUSIONS: Reinforcement learning offers several distinct advantages for community detection, including scalability and knowledge of the walk trajectories defining those communities. Applied to currently available human protein interaction networks, this method had comparable accuracy with other algorithms and notable savings in computational time, and in turn, led to clear predictions of protein function and interactions for several uncharacterized human proteins.

Integrated modeling of the Nexin-dynein regulatory complex reveals its regulatory mechanism.

Cilia are hairlike protrusions that project from the surface of eukaryotic cells and play key roles in cell signaling and motility. Ciliary motility is regulated by the conserved nexin-dynein regulatory complex (N-DRC), which links adjacent doublet microtubules and regulates and coordinates the activity of outer doublet complexes. Despite its critical role in cilia motility, the assembly and molecular basis of the regulatory mechanism are poorly understood. Here, utilizing cryo-electron microscopy in conjunction with biochemical cross-linking and integrative modeling, we localized 12 DRC subunits in the N-DRC structure of Tetrahymena thermophila . We also found that the CCDC96/113 complex is in close contact with the N-DRC. In addition, we revealed that the N-DRC is associated with a network of coiled-coil proteins that most likely mediates N-DRC regulatory activity.

Estimating error rates for single molecule protein sequencing experiments.

The practical application of new single molecule protein sequencing (SMPS) technologies requires accurate estimates of their associated sequencing error rates. Here, we describe the development and application of two distinct parameter estimation methods for analyzing SMPS reads produced by fluorosequencing. A Hidden Markov Model (HMM) based approach, extends whatprot, where we previously used HMMs for SMPS peptide-read matching. This extension offers a principled approach for estimating key parameters for fluorosequencing experiments, including missed amino acid cleavages, dye loss, and peptide detachment. Specifically, we adapted the Baum-Welch algorithm, a standard technique to estimate transition probabilities for an HMM using expectation maximization, but modified here to estimate a small number of parameter values directly rather than estimating every transition probability independently, which should help prevent overfitting. We demonstrate a high degree of accuracy on simulated data, but on experimental datasets, we observed that the model needed to be augmented with an additional error type, N-terminal blocking. This, in combination with data pre-processing, results in reasonable parameterizations of experimental datasets that agree with controlled experimental perturbations. A second independent implementation using a hybrid of DIRECT and Powell's method to reduce the root mean squared error (RMSE) between simulations and the real dataset was also developed. We compare these methods on both simulated and real data, finding that our Baum-Welch based approach outperforms DIRECT and Powell's method by most, but not all, criteria. Although some discrepancies between the results exist, we also find that both approaches provide similar error rate estimates from experimental single molecule fluorosequencing datasets.

Rapid, scalable, combinatorial genome engineering by marker-less enrichment and recombination of genetically engineered loci in yeast.

A major challenge to rationally building multi-gene processes in yeast arises due to the combinatorics of combining all of the individual edits into the same strain. Here, we present a precise and multi-site genome editing approach that combines all edits without selection markers using CRISPR-Cas9. We demonstrate a highly efficient gene drive that selectively eliminates specific loci by integrating CRISPR-Cas9-mediated double-strand break (DSB) generation and homology-directed recombination with yeast sexual assortment. The method enables marker-less enrichment and recombination of genetically engineered loci (MERGE). We show that MERGE converts single heterologous loci to homozygous loci at ∼100% efficiency, independent of chromosomal location. Furthermore, MERGE is equally efficient at converting and combining multiple loci, thus identifying compatible genotypes. Finally, we establish MERGE proficiency by engineering a fungal carotenoid biosynthesis pathway and most of the human α-proteasome core into yeast. Therefore, MERGE lays the foundation for scalable, combinatorial genome editing in yeast.

Amino acid sequence assignment from single molecule peptide sequencing data using a two-stage classifier.

We present a machine learning-based interpretive framework (whatprot) for analyzing single molecule protein sequencing data produced by fluorosequencing, a recently developed proteomics technology that determines sparse amino acid sequences for many individual peptide molecules in a highly parallelized fashion. Whatprot uses Hidden Markov Models (HMMs) to represent the states of each peptide undergoing the various chemical processes during fluorosequencing, and applies these in a Bayesian classifier, in combination with pre-filtering by a k-Nearest Neighbors (kNN) classifier trained on large volumes of simulated fluorosequencing data. We have found that by combining the HMM based Bayesian classifier with the kNN pre-filter, we are able to retain the benefits of both, achieving both tractable runtimes and acceptable precision and recall for identifying peptides and their parent proteins from complex mixtures, outperforming the capabilities of either classifier on its own. Whatprot's hybrid kNN-HMM approach enables the efficient interpretation of fluorosequencing data using a full proteome reference database and should now also enable improved sequencing error rate estimates.

Native doublet microtubules from Tetrahymena thermophila reveal the importance of outer junction proteins.

Cilia are ubiquitous eukaryotic organelles responsible for cellular motility and sensory functions. The ciliary axoneme is a microtubule-based cytoskeleton consisting of two central singlets and nine outer doublet microtubules. Cryo-electron microscopy-based studies have revealed a complex network inside the lumen of both tubules composed of microtubule-inner proteins (MIPs). However, the functions of most MIPs remain unknown. Here, we present single-particle cryo-EM-based analyses of the Tetrahymena thermophila native doublet microtubule and identify 42 MIPs. These data shed light on the evolutionarily conserved and diversified roles of MIPs. In addition, we identified MIPs potentially responsible for the assembly and stability of the doublet outer junction. Knockout of the evolutionarily conserved outer junction component CFAP77 moderately diminishes Tetrahymena swimming speed and beat frequency, indicating the important role of CFAP77 and outer junction stability in cilia beating generation and/or regulation.

Does AlphaFold2 model proteins' intracellular conformations? An experimental test using cross-linking mass spectrometry of endogenous ciliary proteins.

A major goal in structural biology is to understand protein assemblies in their biologically relevant states. Here, we investigate whether AlphaFold2 structure predictions match native protein conformations. We chemically cross-linked proteins in situ within intact Tetrahymena thermophila cilia and native ciliary extracts, identifying 1,225 intramolecular cross-links within the 100 best-sampled proteins, providing a benchmark of distance restraints obeyed by proteins in their native assemblies. The corresponding structure predictions were highly concordant, positioning 86.2% of cross-linked residues within Cɑ-to-Cɑ distances of 30 Å, consistent with the cross-linker length. 43% of proteins showed no violations. Most inconsistencies occurred in low-confidence regions or between domains. Overall, AlphaFold2 predictions with lower predicted aligned error corresponded to more correct native structures. However, we observe cases where rigid body domains are oriented incorrectly, as for ciliary protein BBC118, suggesting that combining structure prediction with experimental information will better reveal biologically relevant conformations.

Label-free proteomic comparison reveals ciliary and non-ciliary phenotypes of IFT-A mutants.

DIFFRAC is a powerful method for systematically comparing proteome content and organization between samples in a high-throughput manner. By subjecting control and experimental protein extracts to native chromatography and quantifying the contents of each fraction using mass spectrometry, it enables the quantitative detection of alterations to protein complexes and abundances. Here, we applied DIFFRAC to investigate the consequences of genetic loss of Ift122, a subunit of the intraflagellar transport-A (IFT-A) protein complex that plays a vital role in the formation and function of cilia and flagella, on the proteome of Tetrahymena thermophila . A single DIFFRAC experiment was sufficient to detect changes in protein behavior that mirrored known effects of IFT-A loss and revealed new biology. We uncovered several novel IFT-A-regulated proteins, which we validated through live imaging in Xenopus multiciliated cells, shedding new light on both the ciliary and non-ciliary functions of IFT-A. Our findings underscore the robustness of DIFFRAC for revealing proteomic changes in response to genetic or biochemical perturbation.

Mining Patents with Large Language Models Elucidates the Chemical Function Landscape.

The fundamental goal of small molecule discovery is to generate chemicals with target functionality. While this often proceeds through structure-based methods, we set out to investigate the practicality of orthogonal methods that leverage the extensive corpus of chemical literature. We hypothesize that a sufficiently large text-derived chemical function dataset would mirror the actual landscape of chemical functionality. Such a landscape would implicitly capture complex physical and biological interactions given that chemical function arises from both a molecule's structure and its interacting partners. To evaluate this hypothesis, we built a Chemical Function (CheF) dataset of patent-derived functional labels. This dataset, comprising 631K molecule-function pairs, was created using an LLM- and embedding-based method to obtain functional labels for approximately 100K molecules from their corresponding 188K unique patents. We carry out a series of analyses demonstrating that the CheF dataset contains a semantically coherent textual representation of the functional landscape congruent with chemical structural relationships, thus approximating the actual chemical function landscape. We then demonstrate that this text-based functional landscape can be leveraged to identify drugs with target functionality using a model able to predict functional profiles from structure alone. We believe that functional label-guided molecular discovery may serve as an orthogonal approach to traditional structure-based methods in the pursuit of designing novel functional molecules.